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Albatrosses can travel a thousand kilometers daily over the oceans. This feat 
is achieved through dynamic soaring, a non-flapping flight strategy where 
propulsive energy is extracted from horizontal wind shears. Dynamic soar¬ 
ing has been described as a sequence of half-turns connecting upwind climbs 
5 and downwind dives through the surface shear layer. We analytically and nu¬ 

merically investigate the aerodynamically optimal flight trajectory for vary¬ 
ing shear thicknesses. Contrary to current thinking, but consistent with GPS 
recordings of flying albatrosses, in thin shears the optimal trajectory is com¬ 
posed of small angle arcs. Essentially, the albatross is a flying sailboat, sequen- 
10 tially acting as sail and keel, and most efficient when remaining crosswind. 

Our analysis constitutes a general framework for dynamic soaring, and more 
broadly energy extraction in complex winds. 
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Dynamic soaring (DS) is the flight technique where a glider, either a bird or man-made, ex¬ 
tracts its propulsive energy from non-uniform horizontal winds such as those found over the 
15 oceans (1). Albatrosses, the archetypal dynamic soarers, have been recorded to travel 5000 km 
per week while relying on wind energy alone (2-4). The potentialities of man-made DS systems 
are humbling: a robotic albatross could patrol the oceans and collect oceanic and atmospheric 
data, traveling at over 40 knots with a virtually infinite range (5, 6). Another flavor of dynamic 
soaring, harvesting the wind shear of the jet stream, is an active field of research (7). 

20 Despite important efforts (5, 8-15), a DS robotic system has remained out of reach. A major 
obstacle to man-made DS has resided in the complexity of the wind power extraction process 
that, by nature, requires planning on the go an energy positive trajectory in a stochastic, hard 
to measure, and poorly understood wind field. Until now, the fundamentals of the energetic 
exchange in DS have only been partially understood, preventing us from efficiently solving 
25 the planning problem. Here, we provide a surprising and yet very intuitive description of DS, 
backed up by analytic and numerical analysis as well as comparison with field data. 

In the first attempt to describe DS, Rayleigh (i) modeled the wind profile as a still boundary 
layer separated from the above windy free stream blowing at Wq by an infinitely thin shear layer 
(see Fig. [^, hereafter Rayleigh’s wind model). He noticed that when traversing the shear layer 
30 directly up- or downwind, the albatross’ groundspeed is conserved but its airspeed is not, and 
may increase by up to Wq. Rayleigh connected up- and downwind transitions with half-turns in 
order to construct an energy neutral trajectory (hereafter Rayleigh’s cycle, Fig.[^ and e.g. (16)): 
at each transition, the airspeed gain compensates the inherent losses due to drag. Because the 
drag is quadratic with airspeed, a limit cycle is reached. This description of the DS trajectory 
35 has carried on until today (5, 10-19) in two energetically equivalent forms: trajectories with 
constant turn direction are 0-shaped, or loitering; trajectories with alternating turn directions 
are S-shaped, or traveling. 
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Figure 1: Wind profile. (A) Wind field behind waves. Color-eoding: wind intensity. Experi¬ 
mental data adapted from (20). (B) The logistie wind profile in this study eaptures adequately 
the wind field in separated regions, sueh as behind oeean waves. More generally, it eonstitutes a 
robust way to approximate a wide elass of wind fields, based on two parameters: a typieal wind 
speed inhomogeneity Wq separated by a buffer of typieal dimension 6 . (C) Rayleigh’s wind 
model is the limit of the logistie profile for 5 —)■ 0. 
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Recent observations based on high-accuracy GPS measurements {2,21) (reproduced in Figs.[^ 
and[^ show that albatrosses do not follow half-turns, but rather an elongated, albeit oscillating, 

40 trajectory. Statistical analysis of this data shows that when flying crosswind albatrosses typi¬ 
cally only turn by about 55°, a mere third the half-turn’s 180°. In this article we explain how 
fundamental this distinction is. We numerically studied the dependence of the aerodynamically 
optimal trajectory of DS on shear layer thickness. We discovered that contrary to prevailing 
theory, in the thin shear layer regime it is a sequence of arcs of vanishingly small angle, with 
45 the direction of flight nearly cross wind at all times. We were able to explain this observation 
analytically, lowering the wind required for DS by over 35% compared to previous models (5). 
For the albatross, the half-turn picture with up- and downwind transitions is misleading, as it 
is suboptimal both energetically and for travel speed. Our theory conceptually unifies dynamic 
soaring, gust soaring (22), turbulence soaring {23, 24), and other wind energy harvesting tech- 
50 niques such as sailing. 

The energy extraction mechanism in DS relies on a transfer of momentum from fast to slow 
air, and any theory starts with formulating the structure of the wind field. In the last two decades, 
a popular approach has consisted in attempting to perform accurate numerical modeling of the 
albatross flight in logarithmic or power law profiles, deemed good models of the average wind 
55 field in the first 20 m above water , where the albatross flies. However, in this framework it has 
been shown {8, 25) that DS is extremely sensitive to the wind field in the first meter above the 
surface, precisely where wind-wave interactions and temporal variability make the logarithmic 
model less accurate. 

In contrast, Rayleigh’s wind model has merit beyond the realm of qualitative analysis for 
60 modeling the sharp wind shear in separated regions, such as behind breaking waves or mountain 
ridges. Recent studies suggest that wind separation in ocean wave fields may be more frequent 
than previously believed ( {20, 26, 27) and Fig. [T]4), further reducing the relative merit of log- 


4 




B C 


Continuous 
airspeed loss 
due to dragy 


Upwind jet 
in wind layer 


/ 

Instantaneous \ 


^ l\ ' 


at Iny-ir transition 




^Continuous 
airspeed loss 
due to drag 

’/'O 

Downwind jet 
in boundary layer 

Top view 

Boundary layer 

Top view 

Boundary layer 



Figure 2: The albatross’ trajectory. (A) Recording of a flying albatross from {21) (top view). 
In crosswind flight the typical turn of the albatross is less than 60° {c.f. ST4). Orange portions 
of the trajectory: the albatross is involved in a 60 ± 20° turn. Red portions: the albatross is 
involved in a 60 ± 10° turn. Note that while in the ground frame the mean albatross travel has 
a downwind component, in the frame moving with the average wind it is nearly crosswind. (B) 
The Rayleigh cycle describes the albatross’ flight as a sequence of half-turns between the windy 
and slow regions. At each layer transition, there is an airspeed gain equal to the wind speed, 
which compensates inherent drag losses that are quadratic in airspeed. However this trajectory 
is suboptimal for energy extraction. Instead, the optimal cycle (C) is composed of a succession 
of small angle arcs. The flight portion in the wind layer is functionally analogous to the sail of 
a sailboat while the portion in the slow layer is analogous to the keel of a sailboat (D). 
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based approaches. 

We modeled the wind with a logistic profile (Fig.[T^) parameterized by the free stream wind 
65 speed Wq and the shear layer thickness 6 

Wo 


W{z) = 


( 1 ) 


1 + exp —z/6 

This captures not only the main features of separated winds over ocean waves but more gener¬ 
ally of any flow with a typical wind inhomogeneity Wq developing over a typical length-scale 
6, such as in turbulence soaring. The regions 2 ; —6, \z\< A5, z ^ 5 represent the boundary 

layer or separated region, shear layer and windy free stream layer, respectively. In the thin shear 
70 limit 5 —)■ 0 the model converges to Rayleigh’s. 

We numerically addressed the following question: for a given shear thickness 5, what is 
the trajectory that requires the least amount of wind? We applied a direct collocation approach 
to a 3D point mass glider model defined by its wing loading m/S and lift-drag coefficient 
curve Cl cd{cl), in particular its glide ratio {cl/cd)^^^^ (c-f- STl and (5, 9, 25)). The 
75 model non-dimensionalization shows the importance of two parameters: the glider’s cruise 


speed Vc = and its associated length A = /g. Typical values for the wandering 
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\pS 

albatross are 14 = 15 m/s and A = 22 m {28). 

Starting from thick shear (5 3> A) and reducing it progressively until the thin shear regime 
{5 A) was reached, we computed the optimal loitering (circling and therefore constrained to 
be half-turn based) and traveling trajectories (Fig. and ST2). 

In thick shear all trajectories are significantly three-dimensional, the loitering and traveling 
trajectories are quantitatively similar, and the turn amplitude of the traveling trajectory is large. 
In thin shear, the loitering and traveling trajectories are qualitatively different. While the loi¬ 
tering trajectory remains significantly 3D, the traveling trajectory’s extension in the 2 ;-direction 
85 shrinks and it becomes approximately 2D. For very thin shears it approaches a quasi-straight 
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Figure 3: Minimum wind trajectories for three shear layer thieknesses {c.f. ST2). On the left, 
the trajeetories are eonstrained to fulfill the specific requirement that the heading increases by 
360° over a cycle, hence their loitering appearance. On the right, the heading is required to be 
periodic, hence their traveling appearance. For the 3D trajectory the scale is common and is 
indicated on the bottom right comer: the trihedral is of length A = ^ (22 m for an albatross). 
Similarly, the scale bars on the top views are of length A. The middle plots 5/\ = 1/128 
are representative of the shear thickness experienced by albatrosses. The traveling trajectory 
requires less wind than the loitering one, with an increasing advantage for thinner shears. When 
5 —)■ 0, the traveling trajectory becomes 2D and is composed of a sequence of vanishingly 
small arcs of finite curvature performed at nearly constant speed. The behavior of the loitering 
trajectory is qualitatively different: for decreasing shear thicknesses, it quickly converges to a 
limit trajectory that remains significantly 3D even for an infinitely thin shear layer. 


7 













line, composed of zigzags of only a few degrees in amplitude. Most importantly, it requires 
only about 2/3 as much wind as the loitering trajectory (Fig. |^. The main characteristics of 
our numerical model are strikingly consistent with albatross flight data, especially given the 
uncertainty associated with the wind field. 

90 The convergence of the traveling trajectory to the neighborhood of z = 0 greatly reduces 
the problem complexity and it was possible to build an analytic model in this limit (Fig. 
and STl). The cycle may be decomposed into glide phases on either side of, but close to, 
z = 0 where the wind shear is weak and airspeed is lost due to drag, and transitions across 
2 ; = 0 of vanishing duration but finite impulse (called “swoops” in e.g. (22)). Denote xfj the air- 

95 relative heading angle, defined to be 0 when the glider is flying crosswind. During glide the 3D 
equations under 2D constraint simplify to dV/d\^p\ = 4 j 2 ^ 4 - During the transition 

through z = 0, for large glide ratios the change in airspeed is AV = fFosin |^o|- Balancing 
the airspeed loss during glides and airspeed gain at transitions brings the relation between wind 
intensity Wq, average airspeed V and heading angle at transition i/jq: 


1 / 


sin-00 Cl 




- V^/c\V^ ‘2cd 

100 The minimum wind-airspeed pair 

W* = 


( 2 ) 


( 3 ) 


105 


TVco'“’ ■ VST 
is reached for -00 —)■ 0 (small amplitude turns), not for 'ipQ = 7r/2 (half-turns). This represents 
the DS analogous to the famous Betz limit of wind energy (29). Furthermore, it can be shown 
that it constitutes the benchmark for soaring in transverse turbulence when the characteristic 
turbulence scale is small compared to A. Lastly, our analytic model is consistent with our 
numerical results (Fig. |^, and consequently provides strong qualitative agreement with the 
albatross trajectory recordings. 
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Figure 4: Minimum wind and turn amplitude of the traveling and loitering trajeetories as a 
funetion of the shear thiekness from our numerieal model {c.f. ST2). (A) For thin shear <5 —>■ 0 
the wind required for the traveling trajectories converges to our 2D model in Eq. Q. Note that 
all glide ratios collapse on nearly the same line. (B) Similarly, the turn amplitude decreases and 
the trajectories become straighten The data from (27) is overlaid: each orange line on the right 
graph represents a turn from the albatross recording in Fig.[^ {c.f. STS and ST4). 
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The previous diseussion was based on airspeed eonsiderations. With the understanding 
gained, we now diseuss DS from an energetie standpoint in the ground frame of referenee 
(Fig. [^). Power is extraeted from the wind through the work done by the lift veetor and is 
110 largest when the glider is flying elose to erosswind. Indeed, the lift veetor is always orthogonal 
to airspeed and flying erosswind allows for a large misalignment between airspeed and ground- 
speed, whieh in turns translates into a large dot produet between the lift veetor and groundspeed. 

The energetie exehange between the air and glider ean be further understood by notieing 
that at eaeh transition, the glider lift ehanges sign. By Kelvin’s eireulation theorem, vortieity 
115 must be shed at eaeh transition. Those vortiees eonstitute the signature of jets formed in the 
wake of the glider whieh manifests a transfer of momentum: the passage of the glider slows 
down the wind layer and aeeelerates the boundary layer. That transfer of momentum is linked 
to an overall reduetion of kinetie energy in the flow as it is harvested by the glider. 

Thus, DS aetually presents strong similarities to sailing: sailboats transfer momentum from 
120 the fast wind to the slow moving oeean by means of two lifting surfaees, the sail and the keel, 
in a manner that is most effleient when operating near erosswind. A DS system effeets a similar 
transfer of momentum from the fast wind to the slow boundary layer, playing sequentially the 
roles of a sail while in the free stream and that of a keel while in the boundary layer (Fig. 
and D). Like the sailboat, it is most effleient when flying near erosswind. 

125 In praetiee, several faetors limit the eyele frequeney of DS. Our analytie model is strietly 
valid for an infinitely large glide ratio in an infinitesimal shear layer. For the albatross, those 
are finite: clIcd < 20 and 5 ~ 1 m. Furthermore, the large albatross wingspan (~ 3 m) ineurs 
a eost of rolling not taken into aeeount in our point mass model, as well as a eonstraint on its 
minimum vertieal travel. Finally, there may be other advantages to lower frequeney eyeles, 
130 sueh as the possibility to synehronize with waves, and a reduetion of the requirement in eontrol 
authority. For the albatross, finite turns are not the eause of energy extraetion, but a eonsequenee 
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of these praetieal limitations. 

This new eoneeptual framework has important implieations: beside eonstituting the basis for 
rieh allometrie extensions, the improved understanding of the albatross flight and the uneovering 
135 of the potentially major role played by wind separation behind waves in its soaring ability help 
refine the eharaeterization and predietion of the albatross’ habitat in a ehanging elimate (30). 
In the quest for a robotie, bioinspired albatross, Eq. Q may well eonstitute the fundamental 
design guideline, while understanding the key roles of shear thiekness and turn amplitude paves 
the way to robust and sealable learning algorithms. 
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STl Analytic model 


Equations of motion and non-dimensionalization 

The analysis utilizes a 3-degree of freedom glider model. Our formulation follows elosely (5,9) 
in the frame or referenee (i, j, k) = (eeast, eNonh, eup) with the state x = (V, ip, 7, z, x, y) where 
5 V is the glider airspeed, -0 is the angle between x and the projeetion of the airspeed V in the 
xOy plane and 7 is the angle between V and the xOy plane and is positive nose up. -0 and 
7 are the air-relative heading angle and air-relative flight path angle respeetively. We assume 
the existenee of a varying wind —W{z)] (blowing from positive to negative y when W > 0). 
The eontrol inputs are the lift eoeffieient and bank angle u = The equations of motion 

10 (eom) are: 


mV = —D — mg sin 7 -|- mW cos 7 sin xp (4a) 

mVp = L cos 0 — mg cos 7 — mbk sin 7 sin ip (4b) 

mV Ip cos j = L sin (p + mW cos ijj (4c) 

i = 1/ sin 7 (4d) 

X = V cos 7 cos Ip (4e) 

y = V cos 7 sin 0 — kk (4f) 


Note that x, y may be considered as output rather than states as they don’t feed back into 
Eqs. ( |4a] - |4d| ). We are looking for the minimum wind trajectories, i.e. periodic in the states 
V, Ip, 7 , 2 ;. In particular we do not require staying upwind. Lift and drag are specified according 
to L,D = 1/2 cl,dpSV‘^. We assume quadratic drag cd = cd,o + kc\ with k~^ = 4/^^Ci5,o 
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15 where /max is the glider maximum lift-to-drag ratio. 

The problem is non-dimensionalized as follows: the natural veloeity of the problem is the 
glider’s eruise speed at cx, = 1, namely K = The natural length is A = {VcY/g. The 

assoeiated time-seale is tc = Vdg = X/Vc. Note that our non-dimensionalization depends only 
on the glider properties and gravity, arguably a more natural ehoiee than approaehes based on 
20 the wind gradient (9). 

Upon non-dimensionalization of the variables v = VjV^ w = Wz = zjX, t = t/tc 
and (■)' = d{-)/dT, Eq. Q beeomes 


v'=— cdv"^ — sinj + w'cosj sind (5a) 

vd = clv^ cos cj) — cos 7 — te' sin 7 sin (5b) 

f COS7/’'= cxn^sin/) + te'cos/' (5c) 

z' = V sin 7 (5d) 

w' = dzwz' (5e) 

x' = V cos 7 cos d (5f) 

y' = V cos 7 sin — tc (5g) 


Two-Dimension limit 


The Rayleigh model is defined as follows: the wind field takes the simplified form w{z) = ^ 
25 for z <0 and w{z) = wq for 2 ; > 0; tco is the difference in wind speed between the two layers. 
z < 0 represents the boundary layer, z > 0 the wind layer, and z = 0 an infinitely thin shear 


layer. The time derivative w' in Eq. (5e) is zero nearly everywhere; the trajectory can therefore 
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be decomposed into an equation of smooth evolution with w' = 0 , hereafter called “glide”, and 
a discontinuous equation of layer transition where w' = wq5 with 5 the Dirac distribution. In 
30 the remaining of this section we drop the subscript, and w = wq is to be understood as the wind 
speed difference between the two layers. 

Assuming that the glider trajectory is also nearly 2-D, the glide and transition equations take 
a simple form, given below. 


Glide 


Consider the dynamics of a glider evolving according to Eq. Q in the vicinity z 
separating plane but not crossing the separation layer. The 2-D approximation z 
7 , 7 ' = 0. Eq. (5b) becomes a constraint on the roll angle cos 0 = and Eq. (j^ 


= 0 ^ of the 
= 0 ^ brings 
simplifies to 


v' = —CdV^ 

( 6 a) 

= CLVsm.(j) 

( 6 b) 


35 Eliminating r, the parametric evolution of v follows: 

dv 1 V 



sign(' 0 ') 


(7) 


reflecting the airspeed cost of turning. The sign function is a consequence of the decrease of 
airspeed with time. Here, / = cl/cd is the glide ratio. 


Layer Transition 

During layer transition the forces remain finite, but w' = w5 induces a finite change of the 
glider’s state. The state transition 1 —>■ can be easily computed from ground- 
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speed eontinuity (a eonsequence of the forees remaining finite). In airspeed quantities it trans¬ 
lates to = V“ ± Wj depending on whether the transition is up or down. This leads to 


tan'ijj^ = tan'll) ± 


w 


V cos'll) 

= v~ \Jl± 2w/v~ sin-^" {w/v~Y 


Note that Eq. (8a) is also smooth near ip = ±7r/2. 


(8a) 

(8b) 


40 Cycle Periodicity 


Both the layer transition and the glide equation are invariant by the transformation (w, 'ip) h-)■ 

{—w, —'ip). This ean be seen as the eonsequenee of the fact that the airspeed gain of flying up¬ 
wind out of the boundary layer is equal to that of flying downwind into the boundary layer. This 
symmetry is particular to the Rayleigh problem: a finite thickness shear layer or a constraint on 
45 the average travel direction would break it. 

As a consequence the physical cycle [transition up—)■ wind layer glide—)- transition down^-boundary 
layer glide] can be subdivided into two equivalent sub-units [transition—)■ glide]—)-[transition—)■ glide], 
expanded below: 


...{'iP^,v^) ^ 

n—1 


glide L 


{'Ip ,v ) 


glide 


{'Ip ,v ) 


n+1 


cycle n 


In a stationary cycle, the airspeed is periodic Vn+i = '^n and the heading angle is anti- 
50 periodic 'ipn+i = —'ipn- Therefore the heading angle evolves by 'ipPj^i —'ippl = ipPj^i -f V'n+i 
a glide phase. 
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Large Glide Ratio Limit 

Previous studies (5) have shown that the necessary wind speed w tends to 0 as the glide ratio (or 
finesse) f tends to oo. We therefore assume / 3> 1 and te -C 1 and look for stationary solution 
55 of the cycle. The transition relation simplifies as follows: to the dominant order '0+ = = ipQ 

and the airspeed gain is Aw = w/v~ sin^o- The glide relation gives the approximated airspeed 
loss during each turn Aw = with a = -^1 — l/c|w^. The anti-periodicity of fio implies 

Afi = 2'0o- Equating airspeed loss and gain brings the equation for the average airspeed 

— , = = —^- w (9) 

Vl - l/ciw4 fio 2 

The minimum w, v pair 

3^/^\/2 q1/4 / /— 

^ — , w = 3 / (10) 

Cl /cd 

60 is attained at the maximum of sine V’o = sinc(O) = 1. Besides, V’o = 0 also maximizes airspeed 
for a given w: the small turn trajectory is optimal both for maintaining airborneness in small 
winds and for maximizing airspeed in large winds. Note that as long as the approximation 
/ S> 1 holds, the glider’s aerodynamic performance measure for the Rayleigh problem is, 
perhaps unsurprisingly, the minimum power coefficient cl^ jcD- 

65 General Case 

Assuming either constant c^, or optimal cl{v) to minimize the airspeed loss, Eq. Q can be 
formally integrated by separation of variables and the exact solution to the Rayleigh problem 
becomes a relatively simple nonlinear algebraic problem. 
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ST2 Numerical solution by direct collocation 


70 


Numerical procedure 


Our numerical model for Fig. j^and|4 is based on the EOM of Eq. j^with 


wiz = 


WQ 


l+exp —zjS' 

We formally rewrite the EOM x = /(x, u). The question that we want answered is the fol¬ 
lowing: For a given glider {cngi /max) a given shear thickness 5, what is the minimum 
wind amplitude Wq that has feasible trajectories, periodic in the state x.^ More specifically, 
75 for the traveling trajectories (right-hand side of Fig. and Fig. Q, the boundary conditions are 
V{T) = 1/(0),'/'(T) = '/>(0),7(T) = 'y{0), z{T) = zff). For the circular trajectories (left-hand 
side of Fig.we imposed the boundary conditions V (T) = V (0),'f(T) = 'f(O) + 2 tt, 'y{T) = 
'y{0),z{T) = z{0),x{T) = x(0). Note that the x-constraint in the latter set of boundary condi¬ 
tions is not strictly required. Without it the upper half cycle tends to peak at a higher altitude, 
80 with very small airspeed and very large cl- The x-constraint maintains to realistic values 
while conserving the main features of the unconstrained trajectories. 

The question is cast into a finite dimensional optimization problem by direct collocation. 
First, time over one period T is discretized into timesteps [0,niT, n 2 T,... ,n 7 v-iT, T] with 
0 <77-1 < < 1- The spacing need not be uniform. We use the shorthand 

85 Xj=x(niT),Uj = u(niT). Following e.g. (25,31), the continuous-time constraints x(niT) = 
u(f))df are approximated by 


Xm, = 7(xi -f Xi_i) - ^(/(Xj, Ui) - /(Xj_i, Ui_i))(ni - ni-i)T 
Z o 


0 = Ci = Xj_i -f - (/(Xi, Uj) -f 4/(xm,, UmJ + /(Xi_i, Ui_i)) (rii - rii^fT 
0 

For the traveling problem, the previous discretization leads to the following nonlinear pro- 
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gram (NLP): 


minimize Wo 

xo,...,xjv,uo,...,ujv,i«o,r 


subject to Cj = 0, i = 1,..., N 


{Vn, ' 4 ’n, 'Jn, zn) — (Vo, 'ipo, 7 o, 2:0) 


( 11 ) 


and zq = 0 


Vi, CL,i > 0 

— 71 < ipi < 71 , — 7 r /2 < 7i < 71/2 

A solution to the NLP is a feasible trajectory that locally minimizes the wind required for flight. 
90 Note that the last three relations are purely technical and the inequalities constraints were not 
active upon solution convergence. 

Similarly, the circular problem is cast into 


minimize wq 

Xo,... ,XJV ,U0,... ,UJV 


subject to Ci = 0, i = l,...,N 


(Vn, 'ipN, 7Af, zn, xn) = (Vo, ^0 + 27 r, 70, zq, xq) 
and Zq = 0 


( 12 ) 


Vi, CL,i > 0 

— Stt < < Stt, —7t/2 < 7* < 7t/2 

The problem was then solved for various (cd,o, /max, 5 ) with a nonlinear solver e.g. SNOPT. 
We typically used N = 140 time steps, leading to 0(1000) variables and constraints. Our 
95 Python implementation converged in 0(1 — 10) minutes on a 2013 Macbook Pro. We used 
more timesteps than in similar studies. The main reason for this choice is that for small 5 the 
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transition through the shear layer is of short duration, and resolving it requires a high level or 
granularity. To reaeh very small values of 6 and validate the eonvergenee of our numerieal model 
to our analytie model, we leveraged on the possibility to utilize non-uniform time spaeing: we 
100 started by solving problems with large 6 and subsequently adressed smaller 6 by adaptively 
refining the time spaeing near the transition in order to maintain a suffieient resolution. 


Results 


The raw results for the eases illustrated in Fig. are eolleeted in Figs. |7]-[^and 10-12 For ease 
with (5 = 1/128 and 1/2048, the eontrol points are non-uniformly spaeed and are more dense 
105 near the transtion ^ = 0. For both the eireular and traveling eases, 6 = 1/128 and 1/2048 are 
qualitatively similar. The boundary thiekness for the albatross is elosest to ease 6 = 1/128. For 
the traveling eases S = 1/128 and 1/2048, the sub-periodieity diseussed in STl is visible - a 
qualitative differenee from S = 2. In eontrast, all eireular eases are qualitatively similar to eaeh 
other aeeross the range of 5’s. 

The main eharaeteristies of the numerieal model are eolleeted in Figs. and 
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ST3 Dimensions for the albatross flight 

Our numerieal proeedure solved the non-dimensional equations of DS. Fig. was dimensional- 
ized with the typieal eharaeteristies of the wandering albatross as used in (28) and reprodueed 
in Table [B 

115 As mentioned in the main text, the exaet strueture of the wind profile is one of the main 
sourees of uneertainty in DS, and the boundary layer thiekness parameter S of our logistie model 
that best represents the wind experieneed the albatross is similarly uneertain. A very reasonable 
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Mass m (kg) 

8.5 

Wing area S (m^) 

0.65 

Glide ratio /max 

20 

Lift coefficient at /max, 

0.5 

X (m) 

21.8 

Vc (m/s) 

14.6 


Table 1: Characteristic of the albatross used in this study. A and Vc are calculated with the air 
density p = 1.2 kg/m^ and acceleration of gravity g = 9.8 m/s^. 

assumption is that the effective shear layer thickness perceived by the albatross must be at least 
of the order of vertical extension of the albatross, from wingtip to wingtip, when it is in a bank. 
120 The albatross’ span being ~3 m, the shear layer perceived by the albatross must be thicker than 
~1 m. Pennycuick’s description of the albatross performing “swoops” at the interface between 
windy and separated regions behind waves suggests that the actual thickness of the shear layer 
is “small”, and that the effective thickness is 0(1 m). Conversely, when the waves are small 
and the wind flow remains attached to the surface, it is possible that the albatross does not 
125 have access to the extremely thin boundary layer and as a consequence perceives a virtually 
thicker shear layer. The published data recording the vertical travel of the albatross is extremely 
scarce. The few data points in (2, 21) suggest that it is in the 5-15m range. Referring to Fig.[^ 
such vertical travels correspond to a shear layer thickness of about 3m and 10 m respectively. 
Therefore, considerations on the albatross size as well as its reported vertical travel suggest to a 
130 shear layer thickness of the order of 1 to 3 meters, and at the very maximum 10 m. In Fig. |^, 
we used the median hypothesis of 2 m with an error bar ranging from 1/2 m to 10 m. 

Note also that the n% shear layer thickness is defined as the thickness, centered around 
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0 , over which the wind field changes by n%, i.e. the height difference Zmin — Zraa.x such that 
W{z^in)/Wo = f % and W{z^^)/Wo = (1 - f)%. 

135 The actual wind during flight of Fig. is reported in (27) to be 7.8 m/s. In Fig. |^, our 
estimate of the wind intensity Wq perceived by the albatross is smaller: as discussed in the 
main text, the albatross is only able to harvest a fraction of that wind difference. Indeed even 
behind separated waves, the mass of air is not at rest with respect to Earth but typically travels 
at e.g. the wave phase speed {32). In non-separated flows the wind at Im is typically more 

140 than 50% the wind at 10 m and here again, the albatross can only exploit a fraction of the total 
wind speed (25). In the present study we assume that the albatross may access 25 to 50% of the 
reported wind speed at 10 m. 


ST4 GPS data analysis 

The data for Figs. and were extracted from the bitmap Fig. 11 of (27) (reproduced in 
145 Fig.[^). For each pixel in the East direction, the center of the trajectory line was determined by 
an average operation. The result was filtered with the filtfilt filter fromscipy.signal. 
The (ground) heading angle was then calculated (Eig.[^). Eigs.[^ and D report the distribution 
of the turns in the recording. The amplitude of each turn (in degrees) is also overlaid in Eig. |^, 
with color darkness proportional to the curvilinear length of the turn (in meters). Similarly, the 
150 beeline progress is overlaid in Eig. |^. 

While our model predicts the albatross’ turn amplitude extremely well, it underpredicts the 
cycle length (Eig. |^). Two factors may explain this: 1) There is a benefit in remaining in the 
“keel” phase of the cycle because the ground effect reduces drag, possible uplift from wave¬ 
generated winds are a secondary source or energy, and the direction of travel is skewed upwind 
155 if more time is spent in the slow layer. This hypothesis is consistent with the observation that 
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Figure 5: Characteristics of the minimum-wind cycle from ST2. Same legend as in Fig. 
(A) Height separation between the lowest and highest point of the eyele. For thin shears the 
traveling trajeetory is nearly 2D. Note that the eonvergenee rate is only about 2; 52 / 3 . (B) 

Cyele duration. (C) Maximum airspeed attained during the eyele. (D) Forward travel during 
one eyele. The orange {resp. eyan) dots correspond to twice the length of the sail {resp. keel) 
phase in Fig. [^. 
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Albatross flight recording 
(Sachs 2016) 



Distribution of the albtross turns 
in (Sachs 2016). Mean turn 55°, std dev 18° 



Turn amplitude (°) 


Figure 6: Analysis of the albatross’ trajectory. (A) Recording of an albatross travelling ac- 
cross a low wind {21). (B) Albatross heading along the trajectory. In (C) the statistical analysis 
of the flight shows that the albatross turns on average by 55° (median also 55°), significantly 
less than the 180° of the half-turn picture. In this particular recording, the albatross virtually 
never turns more than 90°. (D) Curvilinear length of the individual turns. 

the keel phase is indeed typically longer than the sail phase. 2) In conditions of sufficient wind, 
secondary goals of the albatross may include forward travel speed, minimum control activity 
and/or aerodynamic loads and reaching higher altitudes for e.g. a better observation of the 
ocean. Those goals are better satisfied in long cycles. They also tend to skew the trajectories to 
160 smaller turns, consistent with the small over-prediction of the turn amplitude by our model. 
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Figure 7: Solution to the Rayleight problem for /^ax = 20,Ci,/_ = 0.5,5 = A/2. w;o 
0.52. 


26 
















Trajectory top view 





Figure 8: Solution to the Rayleigh problem for /^ax = 20,Ci,/_ = 0.5,5 = A/128. w;o 
0.23. 
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Figure 9: Solution to the Rayleigh problem for /max = 20.ci,/„„ = 0.5, i = A/2048. 

Wq = 0 . 21 . 
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Figure 10: Solution to the Rayleigh problem for /max = 20,cz.j_ = 0.5,5 = A/2. w;o 
0.55. 
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Figure 11: Solution to the Rayleigh problem for /max = 20.ci,/„. = 0.5, i = A/128, 
Wq = 0.304. 
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Figure 12: Solution to the Rayleigh problem for /max = 20, = 0.5,5 = A/2048, 

wq = 0.301. 
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